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Background of Invention 

Field of the Invention 

[0001] The invention relates generally to the field of seismic exploration. More 

specifically, the invention relates to methods for processing seismic data acquired using 
both motion sensors and pressure sensors. 

Background Art 

[0002] In seismic exploration, seismic data are acquired by imparting acoustic energy 

into the earth near its surface, and detecting acoustic energy that is reflected fi"om 
boundaries between different layers of subsurface earth formations. Acoustic energy is 
reflected when there is a difference in acoustic impedance between adjacent layers to a 
boundary. Signals representing the detected acoustic energy are interpreted to infer 
structures and composition of the subsurface earth structures. 

[0003] In marine seismic exploration, a seismic energy source, such as an air gun, or air 

gun array, is typically used to impart the acoustic energy into the earth. The air gun or 
array is actuated at a selected depth in the water typically while the air gun or array is 
towed by a vessel. The same or a different vessel tows one or more seismic sensor 
cables, called "streamers", in the water. Generally the streamer extends behind the vessel 
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along the direction in which the streamer is towed. Typically, a streamer includes a 
plurality of hydrophones disposed on the cable at spaced apart, known positions along the 
cable. Hydrophones, as is known in the art, are sensors that generate an optical or 
electrical signal corresponding to the pressure of the water or the time gradient {dp/dt) of 
pressure in the water. The vessel that tows the one or more streamers typically includes 
recording equipment to make a record, indexed with respect to time, of the signals 
generated by the hydrophones in response to the detected acoustic energy. The record of 
signals is processed, as previously explained, to infer structures of and compositions of 
the earth formations below the locations at which the seismic survey is performed. 

[0004] Marine seismic data often include two particular artifacts that require techniques 

to account for in order to more accurately infer the structure and composition of the 
subsurface earth formations. These two artifacts, known as ghosting and water layer 
multiple reflections, arise because water has a substantially different acoustic impedance 
than the air above the water surface, and because water typically has a substantially 
different acoustic impedance than the earth formations at the bottom of the water (or sea 
floor). 

[0005] Ghosting and water layer multiples can be understood as follows. When the air 

gun or air gun array is actuated, acoustic energy radiates generally downwardly where it 
passes through the sea floor and into the subsurface earth formations. Some of the 
acoustic energy is reflected at subsurface acoustic impedance boundaries between layers 
of the earth formations, as previously explained. Reflected acoustic energy travels 
generally upwardly, and is ultimately detected by the seismic sensors (hydrophones) on 
the one or more streamers. After the reflected energy reaches the streamers, however, it 
continues to travel upwardly until it reaches the water surface. The water surface has 
nearly complete reflectivity (reflection coefficient equal to imity) with respect to the 
upwardly traveling acoustic energy. Therefore, nearly all the upwardly traveling acoustic 
energy will reflect fi-om the water surface, and travel downwardly once again. The water- 
surface reflected acoustic energy will also be shifted in phase by about 180 degrees from 
the upwardly traveling incident acoustic energy. The surface-reflected, downwardly 
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traveling acoustic energy is conunonly known as a "ghosf signal. The ghost signal 
causes a distinct **notch", or attenuation of the energy within a limited frequency, in the 
acoustic energy detected by the hydrophones. The frequency of the notch in the detected 
acoustic signal is related to the selected depth at which the streamer is disposed, as is well 
known in the art. 

[0006] The downwardly traveling acoustic energy reflected from the water surface, as 

well as acoustic energy emanating directly from the seismic energy source may reflect 
from the water bottom and travel upwardly, where it is detected by the hydrophones. 
This same upwardly traveling acoustic energy will also reflect from the water surface, 
once again traveling downwardly. Acoustic energy may thus reflect from both the water 
surface and water bottom a number of time before it is attenuated, resulting in so-called 
water layer reverberations. Such reverberations can have substantial amplitude within the 
total detected acoustic energy, masking the acoustic energy that is reflected form 
subsurface layer boundaries, and thus making it more difficult to infer subsurface 
structures and compositions from seismic data. 

[00071 It is known in the art to provide a so-called "dual sensor" cable for detecting 

acoustic (seismic) signals for certain types of marine seismic surveys. One such cable is 
known as an "ocean bottom cable" (OBC) and includes a plurality of hydrophones 
located at spaced apart positions along the cable, and a plurality of substantially 
collocated geophones on the cable. The geophones are responsive to the velocity of 
motion of the medium to which the geophones are coupled. Typically, for OBCs the 
medium to which the geophones are coupled is the water bottom or sea floor. Using 
signals acquired using dual sensor cables enables particularly useftil forms of seismic 
data processing. Such forms of seismic data processing generally make use of the fact 
that the ghost signal is substantially opposite in phase to the acoustic energy traveling 
upwardly. The opposite phase of the ghost reflection manifests itself in the signals 
measured by hydrophones by having opposite sign or polarity as compared with the 
upwardly traveling acoustic energy. Because a geophone is directionally sensitive 
(whereas a hydrophone is not), the downwardly traveling ghost signal causes another 
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phase reversal in the signal as detected by the geophone, with the result that the upwardly 
traveling wave and the downwardly travehng ghost signal detected by the geophone will 
be in phase. 

[0008] The foregoing relationship between polarities of upgoing and downgoing acoustic 

energy has led to a number of "deghosting" and water layer effect attenuation techniques. 
One such technique is described in U. S. patent no. 4,486,865 issued to Ruehle. Pairs of 
detectors each comprise a geophone and a hydrophone. A filter is appUed to the output 
of at least one of the geophone or hydrophone in each pair so that the frequency content 
of the filtered signal is adjusted. The adjustment to the fi-equency content is such that 
when the filtered signal is combined with the signal firom the other sensor, the ghost 
reflections cancel. 

[0009] U. S. patent no. 5,621,700 issued to Moldovenu also discloses using at least one 

pair of sensors in a method for attenuating ghosts and water layer reverberations. 

[0010] U. S. patent no. 4,935,903 issued to Sanders et al. discloses a method for reducing 
the effects of water later reverberations which includes measuring pressure at vertically 
spaced apart depths, or by measuring pressure and particle motion using sensor pairs. 
The method includes enhancing primary reflection data for use in pre-stack processing by 
adding ghost data. 

[0011] U. S. patent no. 4,979,150 discloses a method for marine seismic exploration in 

which output of substantially collocated hydrophones and geophones are subjected to a 
scale factor. The patent states that collocated hydrophones and geophones can be 
positioned at the sea floor or above the sea floor. 

[0012] The benefits of using dual sensor cables have been well recognized. However, 

techniques known in the art for deghosting and multiple attenuation are typically intended 
for use with OBCs. Methods known in the art for deghosting and water layer multiple 
attenuation known in the art work for OBCs because the ghost and water layer multiple 
energy is typically downgoing at the sea floor, making it relatively simple to discriminate 
the ghosts and water layer multiples from seismic energy reflected firom earth structures 
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beneath the sea floor, which is generally upgoing. It is desirable to be able to deghost 
and attenuate water layer multiples in dual sensor signals acquired using streamer-type 
dual sensor cables towed by a vessel. Using dual sensor streamers is particularly 
desirable because moving a streamer from one location to another is much less time 
consuming and much less difficult than moving an OBC. Further, it is desirable to have a 
method for deghosting and multiple attenuation for a dual sensor streamer which is 
relatively insensitive to the water depth at which the streamer is positioned (towed), is 
relatively insensitive to undulations in the water surface, and for which knowing the 
energy source '^wavelet" (acoustic signature) beforehand is unnecessary. However, 
techniques known in the art for deghosting and water layer multiple attenuation have not 
proven very effective for dual sensor streamer data because at the depth in the water at 
which streamers typically are towed, water layer multiples may include both upgoing and 
downgoing components, making directional discrimination difficult. 

Summary of Invention 

[0013] One aspect of the invention is a method for deghosting and water surface multiple 

reflection attenuation in dual sensor marine seismic data. A method according to this 
aspect of the invention includes decomposing the data from a plurality of source positions 
into upgoing and downgoing wavefield components using a measured parameter related 
to pressure and a measured parameter related to vertical particle motion. A multiple-free 
wavefield is then determined from the decomposed wavefields. In this aspect of the 
invention, the deghosting is performed by transforming the data into the spatial frequency 
domain, and separating upgoing and downgoing wavefield components of the 
transformed data. 

[0014] Another aspect of the invention is a method for deghosting seismic data. The 

seismic data comprise measurements of a vertical component of a parameter related to 
particle motion and a parameter related to pressure. The measurements related to 
pressure and particle motion are made collocated, and are made at a plurality of spaced 
apart positions. The method includes transforming the data into the spatial frequency 
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domain, and separating an upgoing wavefield component of the transformed data. The 
upgoing wavefield component is the deghosted wavefield. 

[0015] A method for seismic exploration according to another aspect of the invention 

includes actuating a seismic energy source in a body of water at a plurality of positions. 
A parameter related to pressure is measured at a plurality of locations at a selected depth 
below the surface of the body of water. A parameter related to a vertical component of 
particle motion is also measured at substantially the same locations as measuring the 
parameter related to pressure. The measurements of the pressure related parameter and 
particle motion parameter acquired at each of the plurality of source positions are then 
decomposed into upgoing and downgoing wavefield components, and a substantially 
multiple-firee wavefield is then determined from the decomposed wavefield components. 

[0016] Other aspects and advantages of the invention will be apparent from the following 
description and the appended claims. 

Brief Description of Drawings 

[0017] Figure 1 shows a seismic vessel towing a source and streamer, to acquire signals 

that can be processed according to one aspect of the invention. 

[0018] Figure 2 shows a coordinate system and notation used to explain signal processing 

methods according to the invention. 

[0019] Figure 3 is a flow chart of an example embodiment of a method according to the 

invention. 

[0020] Figure 4 shows a general purpose programmable computer configured to read a 
computer program according to the invention. 

Detailed Description 

[0021] An example technique for acquiring seismic data that can be used with seismic 

data processing methods according to the invention is shown in Figure 1 . The example in 
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Figure 1 is two-dimensional, meaning that the technique is related to data acquired using 
a single seismic sensor cable ("streamed) or an ocean bottom cable. A two-dimensional 
representation of methods according to the invention is only meant to show the principle 
of the invention and is not intended to limit the scope of the invention. Methods of data 
processing according to the invention can also be used with three-dimensional acquisition 
techniques, in which more than one seismic soxirce and/or laterally spaced apart streamers 
(and/or ocean bottom cables) are used to acquire seismic data representative of a selected 
subsiuiace area. 

[0022] In Figure 1, a seismic vessel 10 includes navigation, seismic source control and 

recording equipment (referred to for convenience hereinafter as the "recording system") 
of any types well known in the art and shown generally at 12. The recording system 12 
causes a seismic source 14 towed in the water 1 1 to actuate at selected times. The source 
14 may be any type well known in the art, including air guns or water guns. In other 
embodiments, there may be more than one seismic energy source. The actual seismic 
energy source configuration used is not intended to limit the scope of the invention. 

[0023] In some embodiments, a dual sensor cable or "streamer" 16 is also towed by the 

vessel 10. The streamer 16 includes substantially collocated pairs of sensors 18 at spaced 
apart positions along the cable 16. Each sensor pair 18 includes a sensor (not shown 
separately) responsive to the pressure in the water 11 or to changes in pressure such as 
change in pressure with respect to time, and a sensor (not shown) responsive to particle 
motion of the water 11. As is well known in the art, the pressure responsive sensor may 
be a hydrophone, and the motion responsive sensor may be an accelerometer, a position 
sensor or a geophone (responsive to particle velocity). The type of each of the sensors 
(not shown) actually used in any acquisition system is not intended to limit the scope of 
the invention. For purposes of the invention, and as will be explained below, it is only 
necessary to be able to determine a vertical component of the particle motion (or 
acceleration or velocity) at each sensor pair 18. Sensor arrangements which enable such 
determination include, for example, a geophone having an orientation measuring device 
associated therewith such that the vertical component of motion may be calculated from 
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the actual measured motion and the relative orientation of the geophone. Other examples 
of a sensor include a gimbaled geophone, such that the geophone sensitive axis is 
oriented substantially vertically at all times. A preferred type of streamer including both 
particle motion and pressure sensors is disclosed in U. S. Patent Application No. 
10/233,266, filed on August 30, 2002, entitled, "Apparatus and Method for 
Multicomponent Marine Geophysical Data Gathering", and assigned to the assignee of 
the present invention. Said U. S. Patent Application No. 10/233,266 is incorporated 
herein for reference. 

[0024] In some embodiments, an ocean bottom cable (OBC) shown generally at 23 can 

be deployed on the water bottom 22. The OBC 23 can be of any type well known in the 
art, and includes pairs of substantially collocated pressure and particle motion sensors, 
shown generally at 19, disposed at spaced apart positions along the OBC 23. Signals 
generated by the sensors 19 are recorded by a recording unit 21 for later retrieval and 
processing. In some embodiments, both streamers and OBCs can be used to record 
signals to be processed according to methods that will be explained below with reference 
to Figures 2-4. 

[0025] Still referring to Figure 1, when the seismic energy source 14 is actuated, acoustic 

energy travels downwardly, at 26. Some of the downwardly traveling energy 26 
penetrates the water bottom 22 and reaches a subsurface layer boimdary 24. Acoustic 
energy is reflected from the layer boundary 24, whereupon the reflected energy travels 
upwardly, as shown generally at 30. The upwardly traveling acoustic energy 30 is 
detected by the sensor pairs 18 on the streamer 16 (or the sensors 19 on the OBC 23 if an 
OBC is used). The upwardly traveling energy 28 reflects from the water surface 20, 
whereupon the energy travels downwardly again, as shown at 32. The water surface 
reflected energy 32 is detected by the sensor pairs 18, (and/or 19) resulting in a ghost 
signal, as explained in the Background section herein. The water surface reflected energy 
32 also may be reflected from the water bottom 22, and becomes upwardly traveling 
energy, shown generally at 29. Some of the energy emanating directly from the source 
14 can also be reflected from the water bottom 22 and becomes part of the energy 
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reflected from the water bottom, as shown at 28. Also as explained in the Background 
section herein, acoustic energy can reflect from the water siu-face (downgoing energy 32) 
and can again reflect from the water bottom (upgoing energy 29) a plurality of times, 
resulting in water-layer multiple reflections. 

[0026] As a result of all the foregoing acoustic energy interactions with the water 1 1 and 

the structures below the water 11, the acoustic energy detected by the sensor pairs 18, 
(and/or 19) referred to as a "total wavefield", includes both upwardly traveling energy 
("upgoing wavefield") and downwardly traveling energy ("downgoing wavefield*'). The 
upgoing and downgoing wavefields include components resulting from subsurface 
reflectors, such as boundary 24, and from water surface and water bottom reflections. 

[0027] Figure 2 shows a coordinate system that will be used in the following explanation 
of seismic signal processing methods according to the invention. The coordinate system 
includes an arbitrarily selected origin, shown at O. Position with respect to the origin O 
of any other point in the coordinate system can be defined by a horizontal coordinate Xj 
and a vertical coordinate X3. For convenience, coordinates of a point can be represented 
by the vector notation x . Position of the source (14 in Figure 1) is represented by the 

S R 

notation x , and positions of each of the sensor pairs is represented by the notation x . 
Collectively, the sensor (receiver) pair positions are denoted by the domain representation 

If. 

[0028] The actual wavefield is denoted as in the frequency 

domain (in which frequency parameter S = jco), hi the foregoing expression for the 

actual wavefield, p denotes the pressure and = {Vj,V3} denotes the particle velocity 

(in the vertical and horizontal directions, respectively). The actual wavefield is generated 

by the seismic energy source (14 in Figure 1) located at position x*^ = {x^ ^xl) . The 

actual wavefield is measiired at a selected depth below the water surface 20, where a 
streamer 16 is used. The selected depth will be the total water depth if an OBC (23 in 
Figure 1) is used. The water surface 20 is not necessarily a horizontal plane. The 
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measured actual (or 'total") wavefield includes energy reflected by the earth's subsurface 
geology and the sea surface 20. The sensor pairs 18 on the streamer 16 (and/or OBC 23) 

are located at positions defined, as previously explained, by X gD . The domain of 
locations of the sensor pairs 18 on the streamer 16 (and/or OBC 23) also do not 
necessarily define a horizontal plane. 

[0029] For purposes of the following analysis, the situation is considered wherein the 

R S 

sensor pair depths are below the source (14 in Figure 1) depth, implying that > . 

It is assumed that the water 1 1 in the domain between the streamer 16 and the sea surface 
20 is substantially homogeneous, and has a density represented by p, and a 
compressibility represented by k . The wave speed, c, is given by the expression 

[0030] It is assumed for purposes of the following analysis that there always exists a 

horizontal observation plane 17, disposed at the depth position = xf^^ between the 

water surface 20 and the depth of the streamer 16 (or OBC 23), and that the observational 
plane 17 is also located below the source depth xs , as shown in Figure 2. 

[0031] Some embodiments of a method according to the invention include transforming 

the data acquired at each sensor pair firom the time-space domain into the spatial 
fi-equency domain. One method for transforming data fi'om the time-space domain to the 
spatial fi-equency domain is the Fourier transform. As a matter of convenience, the 
Fourier transform in the horizontal direction (xy) is used. This spatial Fourier transform 
is defined as: 

[0032] uUsa^,x^,s)= j txp{jsa^x^)uix^,x^,s)dA (i) 
[0033] uix^,x^,s) = y- J exp{-jsa^x^) u{x^,x^,s)dA (2) 
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[0034] In the above expressions, represents the horizontal component of the angular 
slowness vector Also in the above expressions, dA represents the elementary area of 
integration. The total wavefield includes both upgoing and downgoing components, as 
shown in Eq. (3) below: 



[0035] 



(3) 



,obs 



[0036] In the observational plane 17, disposed at depth = x^"^ , the decomposition of 

the total wavefield can be carried out in the spatial Fourier domain where the spectral 
counterparts are given by the expressions: 



[0037] 
[0038] 
[0039] 
[0040] 

[0041] 
[0042] 



and 



in which 



(4) 



(5) 



r= 



1 

^2 



2+^1 



Re(r)>0 



(6) 



It is assumed, as shown above in Eq. (6), that the real component of Y is positive. 
Using the field reciprocity theorem (see, for example, J.T. Fokkema and P.M. van den 
Berg, Seismic Applications of Acoustic Reciprocity, Elsevier, Amsterdam, 1993, Chapter 
10.1) it can be inferred that the upgoing wavefield is uniquely detemiined by the 
measured pressure and the normal component of the particle velocity at the sensor 
domain by the expression: 



[0043] 



-1 
2sT 



\ 



(x^ , s) sp exp{Jsa^Xi - sTx^ ) + 
^(x'' , s) expijsa^x^ - sTx^ ) 



v,dA 



(7) 
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[0044] wherejP(x^,5) and V^v^(x^,5) are the frequency domain measured total 

wavefield quantities, namely, the pressxu'e and the nomial component of the particle 
motion (velocity). In the case where the streamer 16 is substantially horizontal the 

normal component ofthe particle velocity becomes the vertical component, V3. Further, 
in Eq. (7), 5f = {Sf^d^ ) denotes spatial differentiation with respect to . 

[0045] The particle velocity and the pressure ofthe upgoing wavefield are related to each 

other by the expression: 

[0046] P (75^^X3, 5) +rp"'' (75^1,^3, 5) = 0 at ^3 = xf"' (8) 

[0047] In the particular case where the domain of sensor locations, , forms a 

R e 

horizontal plane at -^3 — -^3 C^^jt — ^4,3 ), Eq. (7) can be recogmzed as a spatial Fourier 
transform with the result: 

P"''Usa„s) = 

,0048, ^JSt^\^,^Usa,,.ls)-Tp (Jsa^^.s)] 

[0049] Combining Eqs. (5) and (9), the spatial Fourier transform of the upgoing 

wavefield then becomes: 

10050, exp[sr(x3-X3')] ^_^.^^_^^,^^^_r_ y^^^^. (.0) 

[0051] In the particular case where the streamer 16 is substantially horizontal, the depth 

of the observational plane 17 may be set to the receiver/sensor (streamer) depth, whereby 

^3^^ = , and the pressure of the upgoing wavefield can be determined by the 
expression: 
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[00521 {jsa,,x^,s)\ ^^^^ 

[00531 Similarly, the vertical component of the particle velocity can be obtained by the 

expression: 

I00S4I l-\pr,(jsa„x;.^)-rpUsa„x^,s)] C^' 

2/7 

[00551 Note that {p^^^^^^} represents the upgoing wavefield as reflected from the 

subsurface earth geology (e.g. layer 24 in Figure 1), and is obtained from the total 
measured wavefield quantities. As previously explained, Eq, (10) is the spatial Fourier 
transform of the deghosted wavefield. A simple inverse Fourier transform of Eq. (10) 
provides the deghosted wavefield in the time-space domain. 

[00561 The foregoing analysis shows that if both the pressure (or related parameter such 

as pressure gradient) and the particle velocity (or related parameter such as acceleration) 
are measured, it is possible to obtain the receiver-deghosted wavefield without notches in 
the frequency spectrum. The streamer 16 may be positioned at any selected water depth 
and the results are independent of the nature (geometry) of the sea surface 20. Using the 
power reciprocity theorem, (see J.T. Fokkema and P.M. van den Berg, Seismic 
Applications of Acoustic Reciprocity, Elsevier, Amsterdam, 1993, Chapter 10.2) the 
downgoing wavefield can also be uniquely determined by the measured pressure and 

normal component of the particle velocity at the receiver domain ^ by the expression: 



P'"'^''{jsa,,s)= 



[0057] 1 



X eD L 



(x^ , s) sp exp(y5a, jcf + ^rjc^ ) + 



A A (13) 
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[0058] wherep(x^,5) and V^v^(x^,5^) are the frequency domain transforms of the 

measured field quantities, namely, the pressure and the normal component of the particle 
velocity. 

[0059] The vertical component of the particle velocity and the pressure of the downgoing 

wavefield at the depth of the observation plane 17 are related to each other by the 
expression: 

[00601 pv^''^Usa,,x,,s)'rp''''\jsa,^^^^ =0 at ^3 =xf' (14) 

[0061] In the case that the sensor positions within domain define a substantially 

horizontal plane at the depth of the observation plane 17, wherein -^3 — -^3 (^^ — ^^,3 ), 
the expression for the amplitude P*'*^^" reduces to: 

[0062, exp^x?)^. ^ (.5) 

[0063] Combining Eqs. (4) and (15), the spatial Fourier transform of the downgoing 

pressure wavefields becomes: 

.«0«, exp[-srfe-.?)] ^_^.^^_^^,^^^^^_ ^ .^^__^._^^] (16) 

[0065] In the case where the streamer 16 (and/or OBC 23) is substantially horizontal, the 

observational plane 17 may have its depth selected as the streamer (or OBC) depth, 

obs R. 

whereby ^3 = -^3 , the pressure of the downgoing wavefield can be obtained by the 
expression: 

[0066] P"^'"" ^4^s)^'^\p^3 > ^3^ ^s) + rp {jsa^ , , s)] (i 7) 
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[0067] and the vertical component of the particle velocity can be obtained by the 

expression: 

[0068] vt'^'Usa^.x^.s) = ^[^^3 {jsa,,x^,s) +rp {poc,,x^,s)] (ig) 

[0069] Note that the downgoing wavefield {^''^^'^ , V3''^^" } also includes the incident 

wavefield (the energy directly emanating from the source [14 in Figure 1]) and its source 
ghost. Because the observational (or signal measurement) domain is lower in depth than 
the source position x^, the incident wavefield and source ghost components are 
downgoing with respect to the sensor positions. The pressure of the incident ghost 

wavefield is denoted as p^^^'^ . 

[0070] The downgoing wavefield is written as a superposition of the incident wavefield 

and a scattered wavefield p^^^'^^^ ^g: 

[0071] p'""''' = p'"^''' + p"'-""^" for X, = xf' . (19) 

[0072] Note that the incident ghost wavefield p ' vanishes at the water surface 20. 

[0073] If it is assumed that the water surface 20 is a substantially perfectly reflecting 

plane at ^3 = 0 , it can be observed that: 

[0074] p"'(x,M + p'^^x.As) = 0 (20) 

[0075] and because the incident ghost wavefield p'"^'^ vanishes as well at the water 

surface 20, it can be observed that in view of Eq. (19), 

[0076] p"'Usa,As) + p'''-''''^Usa,,0,s) = 0 (21) 

[0077] or, taking into account the propagation path: 

[0078] p""''^'''^im.X3.s) = -exp(-2srx3)^"^(y5ap.X3,5) (22) 
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[0079] Using Eqs. (1 9) and (22), it can be observed that: 
[0080] p'""'^ {jsa^,x^,s) = ^^^^''(>apX3,5) + exp(-2srx3)p"^(y5a,,X3,5) 
at ^3 =^3^'. (23) 

[0081] It can be observed that the wavefields p"^ and p^^^^ in Eq, (23) are the resuUs 

arrived at from the decomposition procedure previously explained above. 

[0082] Moreover, if it is assumed that the wavefield is generated by a point source 

located at X = X , then the incident wavefield, including its source ghost, is given in the 
spatial Fourier domain by the expression: 

sT 

[0084] Li Eq. (24) W^{s) denotes the source wavelet in the frequency domain. 

Therefore, Eqs. (23) and (24) enable determining the source wavelet from the measured 
pressure and particle velocity. A robust manner for determining the source wavelet is to 
minimize the least-square differences between the right-hand sides of Eqs. (23) and (24), 

for all values of OC^ . 

[0085] Having identified the source wavelet and the downgoing wavefield, a next 

processing step is to remove the multiples related to the water sxirface 20. The 
knowledge of the upgoing and downgoing wavefield, together with the source wavelet, 
can be used in water layer multiple removal procedures known in the art. However, it 
will be shown below that in principle the assumption of a planar water surface and 
knowledge of the source including the source wavelet are unimportant for the water layer 
multiple removal procedure and fiirther processing after it. 

[0086] The explanation of multiple removal methods that do not depend on source 

wavelet knowledge and a planar water surface begins with the propagation invariant at 
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the observation level ^3 — xf'^ that follows from the reciprocity theorem. See, for 
example Section 4.2.2 of J.W. Schoolmeesters, Three-dimensional processing of marine 
seismic data by spectral decomposition, Ph.D. Thesis, Delft University of Technology, 7 
June 2001.. In the two dimensional notation shown in Figure 2, the propagation invariant 

at the observation level ^3 = xf*^ is given by the expression: 



[0087] 



p'^Usoc^,x-^,s)vf{-jsa^,x^,s) - 
sa,eR ip\-jsa^,x^,s)v^{jsai,x^,s) 



In y 



(1A = 0 



(25) 



[0088] in which {p^ , } and {p^ , } are two possible acoustic states, A and B. The 

acoustic states will be further explained below with reference to Eq. (28). Writing each 
of the states, A and ^, as a supposition of upgoing and downgoing wavefields, that is: 



[0089] 
[0090] 



-A ^ -A,up ^ -A,doy>n -B ^ -B,up ^ -B.down 



(26) 



and noting that only oppositely propagating waves contribute, Eq. (25) becomes: 



[0091] 



2jv 



say eR 



p'''^{jsa,,x,,s)vr^"{-jsa,,x,,s)- 



—B^down i 



-—B.down A 



i-jsa^ , x^ , s)v^'"'' {jsui ,X2,s) + 



—A, down A 



(,jsa^,x^,s)v^'"Pirjsa^,Xy,s) - 



p'-"''(-jsa,,x,,s)v,'''""^Usa„x„s) 



—Aydown i 



dA = 0 



(27) 



[0092] Using ttie relations for the upgoing wavefields of the type of Eq. (8) and the 

downgoing wavefields of the type of Eq. (14), the second term can be combined with the 
first term, and the third term can be combined with the fourth term in the left-hand side of 
Eq. (27), thus providing the expression: 



[0093] 



2n 1. 



5ai€R L 



p^-'^Usa,,x,,s)vf'''''^i-jsa,,x,,s)- 
p''"'{-jsa„x,,s)v,'''''''\jsa„x„s) 



dA = 0 



(28) 
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[0094] which is the propagation invariant for upgoing and downgoing wavefields. This 

propagation invariant is the basis for the removal of water surface multiples. 

[0095J In one embodiment of a water surface multiple removal procedure, let the 

previously described State A be the desired multiple-free wavefield and let \j>^ 

denote the pressure and particle motion components of the reflected wavefield. The 
reflected wavefield is the wavefield that would occur in the absence of the water surface 
20. The reflected wavefield is upgoing at the depth of the observational plane 17. The 

incident wavefield that generates the reflected wavefield is denoted as {p"'^, Vj"^} . The 

incident wavefield is downgoing at the depth of the observational plane 17, while the 
reflected wavefield is upgoing, as previously stated. The previously described State B is 
the actual total measured wavefield including water surface multiple reflections. State B 

includes an upgoing wavefield {p^^ ^^^^^ and a downgoing wavefield |^''^^''^v^^^^"| 

Decomposing the upgoing and downgoing wavefield are explained above with respect to 
Eqs. (10) and (15). Substituting States A and B into Eq. (28) provides the expression: 



[00961 2^ 



1^ I 



dA = 0 (29) 



[00971 where X-^ = xf^^ . The propagation invariant of Eq. (29) holds, independently of 

the nature of the sea surface 20. In addition, the propagation invariant holds for any 

choice of incident wavefield that generates a reflected wavefield , A consequence is 
that it becomes possible to choose a source position and a source wavelet arbitrarily. The 

multiple-removal procedure should then lead to the multiple-free wavefield {p'^ ^ 

Therefore, in one embodiment a point source position is selected at x'^'^ with a selected 

desired wavelet W^^{s) , so that: 
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[0099] To make the removal procedure operational, the propagation invariant is 
considered in the spatial domain. Using Parseval's theorem, Eq. (29) becomes: 



[001001 I 

;c,eR 



_p"^{x^,x-^,s)%"'(x^,x^,s) 



dA = 0 (31) 



where = xf^^ . In order to indicate the different source positions X*^ of the actual (measured) 

Sr 

wavefield and X of the desired, multiple-free wavefield, the propagation invariant can be 
written as: 

J p'(x I x*,j)v3*'""(x I x^^)dA = 
1001011 "jV(x|x',.)vr(x|x-,.)dA '■''^ 

[00102] The right-hand side of Eq. (32) is known from the upgoing and downgoing 
wavefield decomposition, as previously explained. In the left-hand side of Eq. (32), 

-rfow/i while p'^is unknown, and represents the multiple-free wavefield. The 

multiple-free wavefield can be determined by solving a system of equations. The system 
of equations in one embodiment includes providing the equivalent of Eq. (32) for each of 

S Sr 
a plurality of source positions X around one particular point source position x . The 

Sr 

particular point source position x may be selected arbitrarily. 

[00103] When the point source is xised as a source position for the desired wavefield, in 
one embodiment, the equations are transformed to the spatial frequency domain by using 
Parseval's theorem in the right-hand side of Eq. (32), substituting Eq. (30), and reverse 
transforming to the space domain. Then, the result is: 
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X, €R 

[00104] 1 (33) 

[00105] Note that, after the calculation of the pressvure p'*(x | X'^'',^) , the related 

particle velocity V3(x|x^'^,5) in the substantially multiple-free reflected wavefield 

follows simply from Eq. (8), by using the previously described spatial Fourier transform 
and inverse Fourier transform. 

[00106] The result of Eq. (33) shows one embodiment of the water surface multiple 
removal procedure. Although Eqs. (28) and (29) are formulated in the spatial and 
frequency domain, these equations can also be written in the spatial Fourier and 
frequency domain. 

[00107] For the assumed case of a horizontal, planar streamer (or OBC), the processing 
sequence in the frequency domain can be simmiarized as follows and as shown in the 
flow chart in Figure 3. The processing begins, at 40, with the measured data of the 
pressure, and at least the vertical component of the particle velocity, 

V3(x^ l^^j-^) where x*^ is the source position for each actuation of the seismic 
source. The source label may also include a position corresponding, for example, the 

midpoint of a source array where an array of air guns is used, x is an element of the 
receiver array. 

[00108] For the actual removal procedure, first the upgoing wavefield, 
{^"^(x^ |x^,^Xv3^(x^ Ix^,^)}, using the spatial Fourier transform and Eqs. (11) 

and (12), and the downgoing wavefield, {^^^^(x^ | x'^,5), vf^^{x^ | X'^,^')} , using 

Eqs. (17) and (18) are decomposed, as shown at 42 . In fact, the upgoing wavefield is the 

deghosted wavefield, as previously explained and as shown at 44. This processing step 

can be carried out for each source label separately (shot-based processing). 
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[00109] The water layer multiple removal in one embodiment includes setting up, at 46, 
and solving, at 48, a system of equations for a plurality of source positions (the discrete 
case of Eq. (33)). Advantageously, the invention provides a method for deghosting and 
water surface multiple removal which is substantially independent of the source wavelet, 
the water surface geometry and the depth of seismic sensors in the water. Therefore, 
corrections for undulation of the water surface are not needed, and it is xmnecessary to 
determine the source wavelet using methods according to the invention. Methods 
according to the invention can also effectively deghost and attenuate water layer 
multiples in seismic data acquired using dual sensor streamers, which has proven 
impractical using methods previously known in the art. 

[00110] The foregoing embodiments of methods according to the various aspects of the 
invention may be performed by a suitably programmed general purpose computer. An 
example of such a computer is shown in Figure 4 having a central processor 50. The 
processor 50 is coupled to a user input device 54 such as a keyboard, and is coupled to a 
display 52 such as a cathode ray tube (CRT) or flat panel Uquid crystal display (LCD). 
A computer program according to this aspect of the invention may reside on any one of a 
number of types of computer readable medium, such as compact disk 62 insertable into a 
CD reader 56, magnetic "floppy" disk 64 insertable into a floppy disk drive 58, or the 
program may reside in a hard drive 60 within or remote from the processor 50. The 
program includes logic operable to cause a programmable computer to perform the data 
processing sequences described above with respect to Figures 1-3. The particular 
embodiment in which a computer program is stored is not meant to limit the scope of the 
invention. 

[00111] While the invention has been described with respect to a limited number of 
embodiments, those skilled in the art, having benefit of this disclosure, will appreciate 
that other embodiments can be devised which do not depart from the scope of the 
invention as disclosed herein. Accordingly, the scope of the invention should be limited 
only by the attached claims. 
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